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INTRODUCTION 


One  could  easily  Justify  the  study  of  the  distribution  of 
quadratic  forms  from  the  standpoint  that  many  of  the  tests  in 
statistics  are  based  on  the  distributions  of  quantities  which  can 
be  thought  of  as  special  cases  of  either  a quadratic  form  or  func- 
tions of  quadratic  forms.  The  applications  are  too  numerous  to 
mention;  however,  we  shell  list  a few  as  illustrations. 

(i)  The  distribution  of  a definite  quadratic 
form  where  the  components  have  a multi- 
variate normal  distribution. 

(ii)  The  problem  of  finding  the  power  func- 
tion of  the  chi-square  statistic,  for 
large  samples,  can  be  reduced  to  that  of 
finding  the  distribution  of  a positive 
definite  quadratic  form  in  non-central 
normal  variates. 

(iii) The  distribution  of  a form  of  the  serial 
correlation  coefficient  can  be  expressed 
In  terms  of  the  distribution  of  a ratio 
of  two  quadratic  forms . See  Anderson  1 

(iv)  Of  special  importance  is  von  Neumann's 

statistic,  ~9j , £" 20 J,  the  ratio  of 


T~.  Numbers  in  square  brackets  refer  to  bibliography. 
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the  mein  square  successive  difference 
to  the  variance,  used  to  test  whether 
observations  are  independent  or  whether 
a trend  exists. 

(v)  Durbin  end  Batson  (~6  J:  use  a similar 
statistic  to  test  the  error  terms  for 
independence  in  least  squares  regression. 

(vi)  Koopnans  [~\2  says,  "Assuming  a normal 
distribution  for  the  random  disturbance,  the 
mathematical  prerequisite  for  an  estimation 
theory  of  stochastic  processes  is  the  study 

of  the  joint  distributions  of  certain  quadratic 
forms  in  normal  variables".  The  problem 
Kocpman3  considers  is  that  ofl  estimating  the 
serial  correlation  in  a stationary  stochastic 
process . 

(vii)  Tc  test  hypotheses  concerning  variance 
components  in  the  analysis  of  variance,  we 
require  the  distribution  of  an  Indefinite 
quadratic  form. 

(viii)  Hotelling  [~9j , shows  how  the  distribution 
of  the  ratio  of  an  Indefinite  quadratic  form 
In  non-central  normal  variates  to  a definite 
quadratic  form  could  he  used  in  the  theory  of 
selecting  variates  for  use  in  prediction. 
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(ix)  McCarthy  £\ shows  how  the  distribution 
of  the  ratio  of  tvo  definite  quadratic  forms 
could  be  used  to  make  an  F test  in  the 
analysis  of  variance  when  the  assumptions 
of  equal  variances  and  independence  of  the 
observations  are  not  met. 

It  would  make  this  report  quite  lengthy  to  discuss  the 
distributions  of  all  these  functions  of  quadratic  forms.  We 
shall  restrict  ourselves  to  the  study  of  a definite  quadratic 
form  in  both  central  and  non-central  independent  normal  variates, 
giving  an  important  application  for  each  distribution.  Then  we 
shall  discuss  two  special  cases  of  an  indefinite  quadratic  form. 
Finally,  we  shall  discuss  a few  inequalities.  We  give  below  a 
slightly  more  detailed  chapter-wise  breakdown. 

In  Chapter  I we  shall  be  concerned  with  the  distribution  of 
a definite  quadratic  form  in  independent  N(0,1)  variates. 

Robbins  /"l6_7»  has  treated  this  problem  but  we  have  carried  it  a 
bit  further.  Robbins  and  Pitman  £ 17_7>  have  given  an  expression 
for  the  distribution  of  a linear  combination  of  chi-square  variates. 
We  feel  that  we  have  Improved  on  this  form.  We  have  derived  an 
expression  which  depends  only  on  the  value  of  the  determinant  of 
the  form  and  on  the  moments  of  the  inverse  quadratic  form.  The 
expression  is  an  alternating  series  which  converges  absolutely 
and  is  such  that  if  we  stop  after  any  even  power  we  have  an  upper 
bound,  and  if  we  stop  after  any  odd  power,  a lower  bound  to  the 


cumulative  distribution  function.  Hotelling  £ 10_7  and  Our  land 
£8J,  have  suggested  the  use  of  Laguerre  polynomials  in  finding 
distributions  of  quadratic  forme.  A brief  account  of  Hotelling's 
method  will  be  given. 

In  Chapter  II  we  have  derived  an  expression  for  the 
distribution  of  a definite  quadratic  form  in  non- central 
independent  normal  variates  which  depends  only  on  the  value  of 
the  determinant  of  the  form  and  on  the  moments  of  the  inverse 
quadratic  form  in  normal  variates  with  imaginary  means.  This 
statement  will  be  made  clearer  later  on.  This  result  enables  us  to 
find  the  power  function  of  the  chi-square  statistic. 

In  Chapter  III  we  have  discussed  the  distribution  of 
the  difference  between  two  independent  chi-squares  having  different 
numbers  of  degrees  of  freedom.  If  the  degrees  of  freedom  are  the 
same,  the  distribution  becomes  the  same  as  the  distribution  of 
the  sample  covariance  in  sampling  from  a normal  population.  Ve 
have  studied  the  properties  of  this  distribution  in  some  detail. 

In  Chapter  17  we  give  some  inequalities  for  the  distribution 
of  a quadratic  form  in  N(0,  1)  variates  and  also  for  the  general 


case. 
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NOTATION 

All  vectors  are  column  vectors  sad  primes  indioate  their 
transposes. 

"p.d.f,"  stands  for  "probability  density  function"; 
"c.d.f,"  stands  for  "cumulative  distribution  function"; 
"r.v,"  stands  for  "random  variable"; 

"q.f."  stands  for  "Quadratic  form"  ; 

"N(u,a)"  stands  for  "a  r.v.  having  a normal  p.d.f.  with 
mean  u and  standard  deviation  o". 


CHAPTER  I 


THE  DISTRIBUTION  OF  A DEFINITE  QUADRATIC  FORM 
IN  INDEPENDENT  CENTRAL  NORMAL  VARIATES 


1.1  The  problem. 

Suppose  we  hare  a q.f.  Qn  * | Y'A  Y in  Y, , . , . ,YR,  where 

the  Yj^  are  independent  N(0,  1)  variates,  and  where  Y'  = (Yp  ...,Yn), 

where  the  dash  denotes  the  transpose  of  the  column  vector  Y. 

Let  Fn(t)  = PrCQ^  t).  Then  the  problem  is  to  find  FQ(t).  It 

is  well  known  that  we  can  make  an  orthogonal  transformation, 

Y * P I,  say,  where  PP1  » P*P  * I,  such  that  | Y'A  Y = | X'P'APX  * 

1 1 n 2 

■g  X 'DaX  = -g  £ a^,  where  Dft  is  a diagonal  matrix  having  the 

elements  a^,...,an  in  its  main  diagonal,  and  where  a^,..,,an  ere 

the  latent  roots  of  the  matrix  A,  Under  Buch  a transformation, 
Xi„..,Xn  remain  independent  N(0,  1).  So  the  problem  is  now  to 

1 ® a 

find  the  Prfe  £ a.XT  <t),  where  we  assume  that  a.  > 0,  i=l,...,n. 
i=l  “ 1 

1.2  The  solution. 

1 n 2 

Theorem  1.1.  Let  Qn  * ■g  £ a^,  where  the  x ^ ere 


independent  N(0,  1),  and  where  a^  > 0,  i«l, 


,n. 


Let  <£« 


1 n 

7 £ a.  a.  . 
*1=1  1 1 


f 


2 


Then, 


(a)  Pn(t) 


t°/2  »<<>* 

(a  H*°  kI  r(|+k+l) 


* ]£ 

where  E(Q^)  is  the  k-th  moment  of 


(b)  The  aerlee  ie  absolutely  convergent  end  therefore  It 
* is  convergent. 

(c)  For  any  two  positive  integers  r and  s end  every  t > 0, 

2s-2  2r-l 

£ d > F (t)  > £ d.  , where 

k=0  * n k=0  * 


tn/2 

71 

ia1...an) 


*<0* 

r(f+k+i) 


Proof. 


Let  dx  = dx. 


• *1  ^ 2 
,dxn  and  £ B_7  = £\  E ftiXii  *.7* 

1=1  ” 


n 


rn(t) 


<=  (2a) 


f-f 

AJ 


exp  - 


1 n 2 
I L xj  dx. 

* 1 1 


We  shall  make  use  of  the  Dirlchlet  Integral; 


k 


n 


J...J  ft 

^ J J*1 


ff1 

v 4ij 


TTrA 

r/r(^)  + 1.7 

i PJ 


where  -oo  < x1  < oo,  and  f[y  Cy  are  all  positive,  and  the 

n xi  pi  . 

range  of  integration  ia  £(— =)  < 1.  See  Edwards  / 7 /.  If  we 

1 Ci 

expand  the  exponential  in  the  integrand,  we  get 


rn(t) 


(Sir) 


n 

2 


/ 1 «,  2.k 

00  ('  5 f 

s ci a* 

/b_7  **>  « 


We  need  to  evaluate  integrals  of  the  following  type: 

i"'/  (L?)*  ax  . 

/s_7  i 1 ' 


If  we  expand  the  Integrand  according  to  the  multinomial  theorem, 
we  get 


V- 


. +i  =k 
n 


k l 


i 


/•••/ 

A.7 


n 

7T 


21 , 


dx4 


k 


We  now  make  use  of  the  Birichlst  Integral  etated  earlier,  where 

Oj  » (|V/2,  Pj  * 2,  jfj  * (21^+1),  getting 

J 


**% 


J-'f  (Z  X?)  d* 

/s_7  i 1 " («1...»n)1/2r(|+k»i) 


z rd^)...  r(in4|) 


1«  +•  • • +1  *k 

l n 


*1  1n* 

1l^**,in*al  ***an 


The  problem  now  le  to  evaluate  this  last  expression.  Recalling 


2sk 


that  If  X1  la  N(0,  l),  then  E(Xp 


r(|) 


^ , we  find  that 


*(<>*  * * /i  f -? 


a'1'! 


21,  21 

T * T n 

k.»  *1  * • ,An 


i « l i i l 

u.„+i=rr,"V  .i  xn 

in  a,  ...a_ 

1 n 


s 


2-kk.»  s — 

V***+1n“k  Ul 


21.  21 

xx  A...E  ^ 


. 1,1 

n“K  1. .* . . .1  a,A...a_ 

1 n 1 n 


r(ixH|)...r(in4) 


* jt  k!  £ ■ ■ ■ i ' _ _ jT — • 

1l+...+1n»k  t . j . a.1, . .a  n 
1 n 1 n 


. n n 
k+j  £ 

>t)  3t 


E(Oj 

(a1...an)1/2  r(|+k+l)  ^ 


tn/2  J°  (-t)k  B((1n) 

p ( t ) » ‘TTp  ^ V 1 n B 

(a^. . .a^1'2  k»0  kt  r(f+k+l) 


This  proves  part  (a)  of  the  theorem. 

To  shov  absolute  convergence  we  note  that  If 


i&i  > ag  > ...  > aQ  > 0,  then 


2 2 

^ + •••  + r><  sr*xi'  “* 

l n n l 


\ 


n 


Mt.  i i **  o v ^ 

zh  s(i  l h*  - ' 80  that 


.n/2 

r“(t)-I^I775 


00  tk  r(|+k) 

k»0  ki  r(|+k+i)akr(|) 


tn/2  t/an 

-sj-j (g  — ~ — - < oo,  for  finite  t. 


(a^.aj1/2  r(|fl) 


This  proves  part  (b)  of  the  theorem. 

The  bounds  ve  obtain  are  based  on  the  fact  that  if  r and  s are 

2a-2  , _vk 
Z 

k=0 


2b-2  , .k 

any  two  poaitive  integers  > 1,  then  for  z > 0,  Z >a“z  > 

* A 


2r-l  ( vk 
r.  III! 


k«0 


k.' 


This  proves  part  (c)  of  the  theorem. 

In  the  case  where  sene  of  the  latent  roots  are  zero,  l.e. 
when  the  form  is  positive  semi-definite  of  rank  r,  say,  we  need 
only  replaoe  n by  r in  the  theorem  and  in  the  proof. 

Remarks. 

(l)  The  moments  of  Q^,  Z(^)  , are  easy  to  obtain  from  the 
cumulants  of  Q^.  The  r-th  cuaulant  of  is  kr(Q^)»  Z a’r. 
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From  Kendall  £l\ J , we  hare  expreseione  for  the  first  ten  momenta 
in  terns  of  the  cumulants. 

(ii)  Let  S_  he  the  sum  of  the  first  r+1  terms  of  the  series 
for  Fn(t).  Then  SQ,  Sg,  S^,  ...  is  a sequence  of  upper  hounds  and 
S^,  S y S^,  ...  ie  a sequence  of  lower  hounds  to  Fn(t),  If  ®2k+2 
is  the  error  committed  hy  stopping  with  s2k+p  then 


E2k+2< 

• S2r- 

g.i.b. 

S2r+1» 

r*0 ,1, . . . ,k;  k=0, 1( 

E2k+1-  1,u,b 

• s2r  " 

g.l.b. 

S2r-1> 

r*l,2, . . . ,k;  k*l,2 

The  values  of 

1 • u • b • 

and  g.l.b.  S2y<+1 

depend  on  the  values  of 

the  latent  roots . We  note  that  Eg^+g  < Sg^-  S2Jc+1  3 the  last 

term  included,  and  E2k+'' < S2k  “ ®2k-l  = the  ^as“  term  included* 

Hence,  the  error  is  less  than  the  lest  tern  included  and  it  is 
positive  if  we  take  an  odd  number  of  terms  and  negative  if  we 
take  an  even  number  of  terms. 

(ili)  The  above  theorem  seems  to  he  in  several  ways  an 
improvement  over  the  method  given  in  Bobbins  £ 1 6J7. 

1.3  An  application:  The  distribution  of  a sum  of  squares  in 

dependent  variates. 

Suppose  that  X^, , . . ,Xn  have  a Joint  multivariate  normal 
distribution  with  zero  means  end  covariance  matrix  equal  to  A’1. 
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1 n 2 

Then  the  problem  le,  what  is  the  distribution  of  7 £ I.?  Row 

d x 1 

tr£ I I'l  < t 7 . .y  S'"J  JSxp  - 3 I'A  I dx. 

(2«)  x,x  tj 


Make  an  orthogonal  transformation,  I * L Y,  sa j where  L L'  ■ L'L  * I. 

n 


Then  I'X  • Y'Y  and  X’A  X ■ Y'L»A  L Y - L a,Y?,  where  a,,.....aw 

^ 1 1 L i E 

the  latent  roots  of  the  matrix  A.  Rov  make  the  trensf  oraat ion 
■ a^Xp  getting 


are 


Pr/3  X'X  < t J 


a 

2 


(2  «)  /-/  Exp 

£\l  a’1*^  t? 


■g  z'z  dz 


* < ^ we  can  make  use  of 

theorem  1.1. 

Remark. 

Combining  the  results  of  theorem  1.1  and  the  above  application 
it  is  easy  to  show  that  we  could  find  the  distribution  of  a definite 
q.f.,  X'A  X,  in  Xp  , Xn  where  Xp  ...,  XQ  have  a multivariate 

normal  distribution  with  cover  lance  matrix  B-1,  and  this  distribution 

involves  as  parameters  the  oharacterlstlc  roots  of  AB"\ 


We  shall  now  state,  without  proof,  an  obvious  corollary 
to  theorem  1.1,  obtained  by  letting  the  first  latent  roots  be 

a^,  the  next  nig  Intent-  roots  be  ag,  etc. 

Corollary  1.1 


1 , o 2 2 

Let  8 » 2^&l  * r?  + "*  + 8'r  )>  Wrlore  th®  are 

”1  r l 

Independent  r.v. 'a  having  a central  chi-square  distribution  with 


^ my  £ „ 1 g 

m.  degroes  of  freedom.  Let  8’  * £(a,.  X +. ..+  a~  ), 

1 j?  d.  JL  2}.n  X*  IP 


r 

M » Z m, , a.,  > 0,  G (t)  = Pr(S  < t),  then 
1 1 " r 


8r(t) 


...O 


1/2 


oo 

E 

k»0 


*(Sj* 

J. 

r(|+t+i) 


# V # 

where  E(Br)  Is  the  k-th  moment  of  Sr.  The  series  is  absolutely 

convergent,  and  for  any  two  positive  integers  s and  3 end  every 
t>  0, 


2s- 2 2J-1 

E dv  > G„(t)  > E d.  , where 

I— A * * 1 A * 


\ 


Remarks. 


(i)  The  mcaents  of  Sr  are  easy  to  obtain  from  the 

cumulante . The  J-th  cumulant  of  S*  ie  k.(S#)  * £ a a” ^ , 

r Jr  2 wii 

We  can  find  the  first  ten  mcaents  in  terms  of  the  cumulante 
from  Kendall  £ 11.7  ♦ 

(ii)  The  above  corollary  gives  a method  which  seems  to 

be  an  improvement  over  the  one  given  in  Robbins  end  Pitman  £ 17  J . 


s method  of 


-V.H  *t-J» 


In  this  section  we  shall  give  a very  brief  account  of  a 

. 12  2 

method  suggested  by  Hotelling  £lo£/.  Let  ^ • -§(a1^+. . ; 

vhere  the  are  independent  N(0,  l)  variates  and  where  a^  > 0, 
n.  Let  g(q)  be  the  p.d.f.  of  Qn  and  let 


-x  m-1 

f(x)  * — fTa) — * X > 0 * where  m * 2‘ 


Then  the  suggested  expansion  for  g(q.)  is  the  Laguerre  series 


11 


e(<i)  • tO i)  l\ 

r-O  r 

where  L^m“1^(q)i  r»0,  1#...,  la  the  sequence  of  Laguerre  polynomials 
satisfying  the  relation 

/» 00 

J f(l)4”'1)(*)L$"'1)(*)  4x  • t1*'1)  f1} 

0 

1;  J • Ojlyteey  fOF  01  > Oe 


The  Laguerre  polynomial  L^Ia"1^(q)  has  the  following  explicit 

s 

representation : 


(a-1) 


(4)  * 


▼■0 


8»0j 1 1 • • • • 


See  Sgego  £ 1 b£f.  It  follows  from  the  orthogonality  condition  that 


r 00 

J «(«>  4«< 


and  so  h#  is  a linear  funotlon  of  the  moments  of  Qn.  The  series 
«(q)  ■ f(q)  (l-rt»1L1(q)+h2L2(q)+...) 


12 

converges  uniformly  over  the  whole  real  axis  and  so  ve  can  Integrate 
tern  by  term. 

Remarks. 

(I)  The  main  drawback  in  using  the  above  expansion  is 
that  no  convenient  bound  is  known  for  the  error  ccmmltted  by 
stopping  after  a certain  number  of  terms. 

(II)  Hotelling  suggests  that  the  above  series  be  used 
to  find  the  p.d.f.  of  the  ratio  of  a definite  quadrat io  form  to 
a sum  of  squares  and  the  p.d.f.  of  an  indefinite  form  by 
convolution,  since  an  indefinite  form  is  the  difference  of  two 


definite  forms. 


CHAPTER  II 


THE  DISTRIBUTION  OF  A DEFINITE  QUADRATIC  FORM  IN  INDEPENDENT  NON- 
CENTRAL NORMAL  VARIATES 

2.1  The  problem. 

Suppose  we  have  a q.  f,  Y’^Y  in  Yp  Yn  where  the 

are  independent  N^,  1)  variates,  and  where  y*  " (Yp  •••>  Yn). 
Let 

0^%*  •••>  Zj  " Gn(tj 

then  the  problem  is  to  find  Gn(tj  £).  ^t  us  make  an  orthogonal 
transformation  Y ■ TX,  say,  where  PT  ■ TP  • I,  such  that 

n 2 

Y’AY  - Hr 'AFX  - X»DaX  « 2 a^  , 


where  Da  is  a diagonal  matrix  having  the  elements  a^,  aR  In  the 
main  diagonal,  where  a^,  ...,  an  are  the  latent  roots  of  the  mdtrix  A 

If  EY  - £,  then  EX'  * P£  - ^ (say),  where  p - (?p...,5n)  and 

" (p.^,  ♦.«,  M-n) • Hance,  under  suoh  a transformation  X^,  ...,  X^ 

remain  independent  with  the  same  variances  as  Y^,  Yn,  but  the 


mean  of  ^ is  now  ♦ *•#  + Ynl£n)>  where  (r1;L,  Ynl) 

are  the  elements  of  the  i-th  row  of  T * • 

So  the  problem  is  now  to  find  the 

the  are  independent  N(|i^,  1), 

2, 2 The  solution, 

, where  the  X^  are  independent  nCm.^,1) 


Theorem  2,1  Let 


^-1 


n 

Za 


1 ft  „i  g 

and  where  a^  > > . . . > a^  > 0.  Let  Z an  Yfc  , where  the 

k*l 


Y^  are  independent  N(i^^,  l),  i « /T[  , Then 


,-»■  tn/2 


oc 

Z 


U) 


(a)  G (tj  J£)  * — T7?"  2 , 

n (ai***an)1/2  8-0  al  r(|  + s + 1) 


where  og  - E(Q^)8  , and  X ■ •!  Z * 

(b)  The  series  is  absolutely  convergent  and  therefore  it  is  conver- 
gent, 

(c)  For  any  two  positive  integers  r and  k and  every  t > 0, 
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2r 

2 ds  > Gn(tjii)  > 
s-0 


2k-l 

2 d . where 
s-0  8 


o. 


e-x  W2 

(a^» • *an) 


(-t)S 

sj 


r(£  + s + l) 


Proof*  We  know  that  if  the  are  independent  N(p, 1 ) and  if 


1 2 2 


1 n 2 


" f 2 ^i  ’ Y * g 2 Xi  » then  P»d*^»  of  Y is 


n 


-X  -7  ? -1  2° 

e e J y 2 

rrrO 


(Xy 


mi  r(m  + 2) 


Let  /\.7  - £ 2 < t J , d£  - dy±  dyn,  and  X^  « | ^ 


Gn(t}  n) 


s-s 

r*j 


7C  e’Xj  s’*3  y" 1 2° 

5-i 


“*4 


3 h-°  V r(vl: 


then 


dy3 
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..  oo 
e"X  S 


i,  i 

\ 1 \ ** 

\ n 


k-0  V,,+in"k  V-*Vr(V 2),”r(inT 


. I(tjajn)  , 


where 


i(tjajn) 


s~s 

r*j 


n 


n 


ExP  - 2y,  7T  y, 
1 .1-1  0 


i,  -r 


d7j  . 


Expanding  the  exponential  and  making  use  of  the  Diriohlet  integral 
stated  in  (l,2),  we  obtain* 


Gn(t;*i) 


, oo 
“A.  - 

O 2j 


i-i  i 

*1  •••  V. 


k-o  i1+,..+in-k  ...  ini  rci^  |)  ...  r(i  + i) 


n 


oo  / , \r.?+k*w 

2 2 lli;/  w — 

r-0  J1+...+dn"T  3 l*°**^nJ 


rcya^)  ...r(in*jn*~) 

W \ i +j+i 
alX  1 2 •••  an  2 r(k+r+|*l) 
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This  can  be  rewritten  as 


-X  +n/2 


00 

2 c. 


si  r(|  ♦ S + 1) 


, where 


c « 2 2 2 

k-0  i^  + ...+  ifi  * k 


where  d » 


81  r(V  ^(-1>k 

“l1-^  V -V  (Jl-V  r<Vl>-r(Vl> 


The  problem  now  is  to  evaluate 


Let  » * 2 e£ , where  the  Yk  are  independent  N(l^k,  l),  i«  /CI# 


w 

- (2n)’l/2  J Y2rExp  - ^(Y-4t)2  dY 


then 


w 
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Z (?!)  (ii2)^  2r“J  r(r-j+  h 


d-o 


‘2d 


r r(r+  i) 

d-o  3 r(j+  i) 


fc  <-^>3 


5r-d 


, and 


E(Q**)S  - 2’3 


n 


si 


j^+».»+dn"s 


23l  23, 

E T.  X 1 
1 n 

di  d 


n 


n 


2“8  . Z 


di+*»«+dn“s 


81  2S"k 

3^  3n 

dii  • ••*„ 


dl 


“n 


2 in  2 

i -0  i -0 
l n 


h , 


where  h « 


3,  3 r(31*|)...r(3n*|)(^)  ^..tUn) 

*i  *n  rfi^i)  ...  r(ij* |j 


2 \ 


8 

z s 

k-0  1!+..,+!^ 


d 


di+...+dn“s 


i 
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si  xj1...  x>i)k  r(jx+  \)  ...  r(jn+  i) 

...  - - i i * 

ixJ...ini( T(lj+  — )«..r(in+  2^al  ***an 


where  d 


But  this  is  just  cg.  Hence, 


e“\W2 


oo 


G (tjji)  r "TT/P  ^ 

~ (a^.TsTr^2  s«0 


(-t)s 

si 


E(Q^)S 
r(f  ♦ s ♦ l) 


This  proves  (a)  of  the  theorem* 

If  Xr(Qn)  is  the  r-th  cumulant  of  Qr,  it  is  not  difficult  to 
ehow  that 


kill  S 

2 J-l 


Hence,  to  find  the  r-th  cumulant  of 


Qh  we  must  replace  by  ip.^  and 


aj  by  aj1  * Betting,  kr(Q^*)  " 

fore,  c is  the  s-th  moment  of  a 

s 


(r-l)i  ; -r 


I a^(l-rn^)  . 


j=l 


y 


r,v,  whose  s-th  cumulant  is 


There- 
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^ 1 ^ —1  2 

If  Qr  * 2 ^ ai  , where  the  are  independent  N(0,  l). 


then  kr(Q^)  ° 2 a^r  , Therefore,  k?(Q^)  < kr(Q^) 


and  it  follows  that  e(Q^*)s  < E((£)s  < 


r(f  ♦ s) 

W 


Consequently, 


-X  .n/2  oo  .s  r(£  * s) 

On(t!  h)  < S-4-. — - S * 2 


(a1...*n)X/2  -O  81  a8r(f  * 8*1) 


-X  n/2  t/,an 

< — - — — ■ ■—»-  ...  < oo  for  finite  t , 


(a,..,ajl/2  r(|  + 1) 


I n' 


This  proves  (b)  of  the  theorem.  The  bounds  stated  in  theorem  2,1 
are  based  on  the  fact  that  if  r and  s are  any  two  positive  integers 
> 1,  then  for  z > 0, 


i=5ik  > e"z  > 
ki 


2r-l 

2 

k-o 


This  proves  (c)  of  the  theorem. 
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Theorem  2.2,  Uhder  the  conditions  of  theorem  2.1 , 


°n^t;  ®”X  °n(t;  -)  • 


where  Gn(t;  0)  equals  the  Fn(t)  of  theorem  1.1. 


Proof. 

X 

f f n -yA  -1/2  oo  (X  y ) 

Gn(t^  = TT  e 1 e 1 y E dy  . 

n l*J  i-1  1 V0  Vr(Ji+i} 


2Jj+l 


r r 2 

£^-j  |3i  dyi  * Pr  ^"i^alX2J1+l+,‘  ,+anX2Jn+l)^t-7» 


r(-f-) 


where  X*  . is  a.  r.v.  having  a central  chi-equare  distribution 


with  2^+1  degrees  of  freedom.  Hence, 


^1  ^n 

X,  ...x" 

X Tl 


Oa(*iM)  * • 2 E ;~ 

k*0  J,+...+J  «k  "1* ' * * *T.J 
l n 


Pr  /|( aL  X g ji+1+*  • • x 2 ^+1*  < *-7. > 
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and  so, 

On(t|  - e’X  0n^tj  2)  • 

We  note  that  equality  is  attained  for  any  given  t,  by  putting.  \ **  0 
(in  which  case  jx  « 0)  , 

Theorem  2,3.  Crider  the  conditions  of  theorem  2*1, 


Gn(tj  ij)  < On(tj  0) 


Proof. 

Us  inf  a special  case  of  a more  general  unpublished  result  due  to 
Sigeity  Moriguti,  we  find  that 


J*  exp  - i(x-n)2  dx  < J exp  - i x2  dx  . 

fax2<  tj  ~ /'ax2<  t J 2 


We  can  generalize  the  above  inequality  to  the  case  n ■ 2 as  follows t 


u exp  - | ^"(x^^2  ♦ (x2ni2)2_7  dx1  dXg 


2 2 
[ a1x1+a2x2« 


S J exp  - |(x2-u2)2  dx. 


£ *2*2$  C alXl-  t"a2x2-7 


V 
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J s 


tj  fafos  t-&2xzJ 


1 2 

exp  - -x, dx.  exp 
2 _ 2 - 2 1 1 


iW* 


< •/  / 

A ft  V 


£ V2  s *-7/  Vi  * t " a2*2  -7 


12  12 
exp  - exp  -gXgdXg 


/ / 

l a1xl+a2x2  S *-7 


1/2  2»  . . 
exp  - i(x1+x2)  dxL  dx2 


Generalizing  the  method  used  for  the  caee  n=2,  It  is  not 
difficult  to  ehov  that  the  ah ore  inequality  could  be  established 
for  n*3;]tj<>*  • This  completes  the  proof.  Combining  the 

results  of  theorems  2.2  and  2.3  we  have  the  very  useful  inequality: 

Gn (t;0)  % On(tjji)  > e“X  Qn(t;0). 

\ 

The  proof  of  theorem  2.3  was  suggested  to  the  writer  by 
Professor  S.  H.  Roy. 

Remarks. 

(1)  The  introduction  of  an  imaginary  mean  in  theorem  2.1 
is  merely  a mathematical  convenience;  we  could  have  emitted  this 

formulation  altogether  and  merely  stated  that  o_  is  the  s-th 

0 


V 


2k 


moment  of  ft  r.v.  whose  s-th  cumulant  is  ks(Q^*)* 

(ll)  The  moments  c8  are  easy  to  compute  vhen  ve  know  the 
cumulants . Kendall  £ 11_7j  gives  the  first  ten  moments  in  terns 
of  the  cumulants. 

(iii)  It  can  be  shown  that  is  invariant  under 

an  orthogonal  transformation  Y«ri,  say,  such  that  r'A  r * D 

0» 

and  r*r  » r r1  » i. 

2.3  An  application;  The  power  function  of  the  chi- square  statistic. 

Suppose  that  the  observations  from  a random  experiment 
can  fall  into  any  one  of  k oells  and  that  the  expected  number 
of  observations  in  the  i-th  cell  under  and  is  m^  and  m^, 
respectively. 

AEjsife  is.  Hjj!  m^,. . . ,mJc,  and  H^<  m^,...,m^.  Let 


X 


2 

0 


and 


X 


2 


k 

£ 

i-1 


where 


nA  is  the  observed  nuaber  of  observations  in  the  i-th  cell.  If 
1-0  is  the  power  of  the  test,  then 
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0 - ?rf  XgSt  / a^, . . . jin^  J 


. (2,)-  I ■ • S «P  /-  | X2J  fr  4(^i ) 

f *0  £ i7  1-1  *1 


approximately,  for  sufficiently  large  n. . Putting  xA  «*  — , 

vm. 

0 

0 ai"°l 

aj/2  * |j.  * — , and  letting  we  get 

•J  mi 

Pri’x  o as  11  / J “ J oxp/- -|s(y1-iii)f7TTdyi 

ik  n 1 1 

r^v;<i7 

- Ok(t;ji)  . 


Hence,  the  application  of  theorem  2*1  will  give  us  the  power 
function  of  the  chi-square  statistic. 


CHAPTER  III 


THE  DISTRIBUTION  OF  AN  INDEFINITE  QUADRATIC 
FORM  IN  INDEPENDENT  CENTRAL  NORMAL  VARIATES 


3.1  Special  case  1, 

Ae  a preliminary  step  to  the  study  of  the  distribution  of 
an  Indefinite  q.f.,  this  section  Is  concerned  vlth  the  distribution 
of  the  difference  between  two  Independent  chi-squares  having 
different  numbers  of  degrees  of  freedcm.  If  the  degrees  of 
freedom  are  the  same,  the  distribution  becomes  the  same  as 
that  of  the  sample  covariance  in  sampling  from  a normal  population. 
We  shall  study  some  of  the  properties  of  this  distribution  In 
order  to  anticipate  the  behavior  of  the  distribution  of  a more 
general  Indefinite  q.f.  Others  £ 15_7*  £\J»  have 

considered  this  problem  from  a slightly  different  standpoint. 

The  main  results  of  this  section  ere: 

(i)  Recurrence  relations  3,1. 4,  3.1.5*  3.1.6,  3.1.7*  3.1.20 
(11)  Inequalities  3.1.9*  3.1.21 

(ill)  Further  properties  3.1.12,  3.1.13*  3.1.15*  3.1.16 
(iv)  An  application  of  a result  due  to  Berry  £$£},  3*1.19  . 

Ir  Tn.m  * V V *h,r*  *n  “d  \ m Independently 
distributed  with  p.d.f.  hQ(x)  and  hjjx)  respectively,  where 


hnW 


1 


» 


x >0* 


\ 


then  if  the  p.d.f.  of  T njtt  i0 


t>  0, 


t < 0. 


...  that  gmn  the  p.d.f.  *»  *>».  t0  •*»  4t  f0r  ‘<0’  ”Pl“’ 
t ty  -1  end  Interchange  n end  «.  Therefore,  v.  ehall  conelder  only 

the  case  t>  0,  for  definiteness. 

Hence, 


emt  r°° 

5,1,1  gn*"(t>  ’ r(S)r(f)  { 


5-i  5-i 

>2x(x+t)2  x dx,  t>0. 


The  moment  generating  function  of  Tn#tt  10  ~ Se 

n _ 2 

(1-0)"  2 (1*)  2 • M<8>  w Me  that 

(i)  If  n,m  — > oo  bo  that  J — * 1,  then  Tft>a  i« 

asymptotically  normal  (Sjp  *&£)  • 

(il)  if  n — > oo,  tut  a remains  finite,  then  Tn>>  i» 
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asymptotically  normal  (^p) , */n). 

(Ill)  If  m — ► oo.  but  n remains  finite,  then  T„  _ ie 

9 9 n,m 

asymptotically  normal  , */"m) . 

In  3.1.1  let  z/t  » y,  getting 


3.1.2 


8^: 


.<*> 


n-fai 


- 1 


r(f)r(f)  “ 


f 


-2yt 


(i«j )‘ 


?-i 


m 

r 


dy. 


Let  n/2  * p,  m/2  « <i,  and  let 


oo 

3.1.3  I(P,4)  - J e’^U  +y)p  y4  dy. 

0 


Integrate  3.1.3  by  parts  three  separate  times  as  follows: 


(i) 

u - (l+y)p  y4. 

dr  • e~2*^  dy, 

(ii) 

u » (l+y)p  e-2yt. 

dy  - y4  dy. 

(ill) 

u « (l+y)p  y4*1, 

dr  » e"2yV  dy,  getting 

(O'  I(P,1)  - A I(P,1-D  + P I(p-l,q).7, 


m>0,  q>0, 
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(ii)'  I(p,q)  - -(^-)  I(p-l,q+l)  + (~j)  I(P,4+1), 

a >-2,  q+l>0, 

(ill) 1 l(p,q)  * (A)  jf  (q-D  Kp,q-2)  + p i(p-i,»-D 
kt* 

+ 2t(q-l)  l(p,q-l)  + 2tp  l(p-l,q)_7> 
m>2,  q-l>0,  respectively. 

The  restrictions  placed  on  m and  q are  to  prevent  the  Integrated 
part  from  becoming  Infinite.  Equate  (i)*  to  (ill)’  getting 

(iv)  I(p+l,q)  * (|^)  I(p+l,q-l)  + (|^)  I(p,q). 

In  (i)'  replace  p hy  p+1  and  q by  q+1  and  then  substitute  from 
(ii)'  and  (iv)  into  (i)'  getting 

(v)  HS  l(p+l>q+i)  « 2(p+l)(q+l)  I(p,q)  + q(q+l)  I(p+l,q-l) 

+ p(p+l)  I(p-l,q+l). 

Now  q(q+l)  I(p+l,q-l)  + p(p+l)  I(p-l,q+l)  » q(q+l)  I(p-l,q-l) 

+ J/"p(p+l)+q(q+l)-7  I(p-l^+l)  + 2q(q+l)  I(p~l,q),  end  I(p-l,q+l) 
+I(p-l,q)  ■ I(p,q).  Substitute  into  (v)  getting 


V 
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Mr  I(p+l,q+l)  ■ /’2(p+l)(a+l)+p(p+l)+q(q+l)-7  I(P><i) 

♦ q(q.^i)  i(p-i,q-D  +/q(q+i)-p(p+i)_7  Kp-M)* 


so  that  finally, 


'fWJiSftftJ'  Wt)!  t>0- 

The  same  relation  holds  for  t <0  If  we  replaoe  t by  -t  and 
Interchange  n and  m.  Note  that  if  n*m,  the  last  term  rani  shea 
and  the  relation  reduoes  to 

2 

*0+4, n+4^  * (n+l*  ®n*e,n-£(t)  * n(n *2)  ga,u^t 


From  the  first  integration  by  parts  ve  get  the  simple  relation 
5,1,5  Sn-^m+a^  “ & + for  411  *• 

We  now  make  use  of  the  p.d.f, ?s  to  obtain  the  c.d.f. 's. 
For  x>0, 

r00 

■ J « 


s 
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r<|)  r(f)  - S 


oo  oo  n . a . 


oo  n . m . 

-t  r _g—  g"-*- 

Integrate  by  parte  where  dv  * e"  dt,  and  u«  J e J{y+t)  y dy, 


getting 


Pr/ \m  >%J  * 8n,»(l)  + Pr7'Tn-2,m  >J-7' 


Hence,  If  n Is  even. 


(n-2)/2 

3.1.6  fc-/Tn(B>x_7-  £ 'h8re 


e"x 

e2,n(x)  • Br/TSjB>xJ  - . 


If  n Is  odd 


3.1,7  Pr/Tn^>x_7 


(n-3)/2 

«n-2J;m(x)  + Pr^Tl,a>x-7 


Frcn  3.1.6  and  3.1.7  it  ie  eeen  that  if  we  have  a table  of  g^Ct) 
and  if  we  know  Pr/’T^m>  xj  for  all  a,  we  can  find  ?r£ Tn^>  xj 
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for  all  n and  m. 

Our  first  objective,  then,  is  to  be  able  to  find 
8q  m(t)  for  any  n and  m.  If  ve  consider  a rectangular  table 

of  gQ  ffi(t)  having  n columns  and  u rows,  ve  will  find  that 

given  g^  ffi(t)  for  n,m»l,2,3,4,5  and  gg  &(t)  for  m**l,2,.,.,  ve 

can  complete  the  table  using  3.1.4  and  3.1.3* 

First  of  all  ve  can  fill  in  the  second  column  using  the 
fact  that 

e't 

8g#a(^)  * 2^72  » . . . , and 

- 1 

* "7572  • m»l,2,3,4,5  . 

If  ve  let  n=m,  ve  obtain,  as  ve  shall  see  later,  that 
(tj(n-l)/2 

sn,n(t)  * J/2^1)  E(n-l)/2(t)’ 

Letting  n-1,3,5  we  find  that  i K^ft), 

t2 

and  g^  ^(t)«^Kg(t),  vhere  K^ft)  is  the  modified  Bessel  function 

of  the  second  kind  of  order  n.  See  McLaohlan  £ 14_7>  sad 
Watson  £ 21_7. 
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We  shall  now  find  the  p.d.f.  for  cases  where  n«m+2.  If 
ve  complete  the  square  in  3-1.1  end  let  n*m+2,  ve  get 


,t»m  oo  m , oo 

*r*.f  •‘‘V-U  nu 


Using  the  fact  that 


«1/8(|)°  r°°  2 

K (t)  -v — J e ty(y-l)  dy,  and 

^ T(n^)  J 


2'  1 


to  get  that 


m+1 

(|)  2 

^"E(m-l)/2(t)  + K(m+l)/2(t)-7 


If  ve  now  put  m«l,3  in  the  above  expression  ve  find  that 


«3>1(t)  -^/^(t)  + Kx(t)__7,  and 


rolB 


34 


*2 

e5;5(t)  - fK^t)  + Kg(t)_7  . 

We  can  now  use  3.1.5  to  find  g^  ^(t),  g^  ^(t),  g^  ^(t),  and 

gj  ,-(t) . Thera  ere  still  six  values  regaining  to  he  found, 

namely,  g^2(t),  gx^(t),  gj^(t),  g^2(t)  and  g^4(t) 

which  can  he  expressed  In  terns  of  the  incomplete  gamma  function. 
In  fact  we  find  that 


T^W) 


oo 

S' 


n 1 
T-  1 

yy2  dy, 


2'  2t 


and 


t 

e 


r(f)2: 


n , 

r1 


oo  n oo  n 

£ J f S’?  ijJ. 

2t  2t 


Letting  n»3,5  we  get  g^2(t),  g5#2(t),  gj^(t),  and  65jl((t). 
Using  3.1*5  we  obtain  gj^  2(t)  and  g^  ^(t).  Now  we  have  all 
gn  a(t)  for  n,n«l,2,3,4,5,  and  together  with  gg  a(t),  3.1.4  and 
3.1,5  we  can  get  all  the  remaining  m(t). 

Our  second  objective,  then,  is  to  find  Pr/*T,  > x 7 for 

i,m  “ 

all  m.  We  shall  first  give  a method  of  evaluating  this  when  m 


is  even,  and  then  a method  for  any  m. 
Evaluation  of  ?r£ m>  xjf,  when  m Is  even 


Pr^Tl,a>  x-7 


« 1 

“ r(|)r<§) 


00 

/ e-t  H(t) 
z 


dt,  where 


oo 

H(t)  « J*  e'2y 
0 


n i 

■3-1 

j (y+t) 


dy. 


Expend  H(t)  in  a,  Maclaurin  eariee  about  t=0.  Now 

1 00  ^ a.  1\ 

H(k)(t)  - (-l)kn  2r(k^)  f e"2yy2  (y+t)  ^ dy, 

0 


and  H(k)(0)  - (-l)k  k 2 r(k+|)  r(|-k-|)  2 ^ ? 


We  may  write 

r tk  fkl 

H(t)  - L ~ rK;(0)  + F (t),  where 
k-0  **  r 


3 6 


t 

Vt}  * 77  S CU-y)'  I(r+1)  (y)  dy. 

0 


»ot.  ttot  | H(r+1>CT)|  < |H(rtl>(0)|  , iottat 

lKr(t)l  *77377  !h(i,+1)(o)|  . 


Therefore , 


la(r^<o>l 

r(|)r(|)(r+i); 


/>00 

y e*4  tr+1  at. 


Hence,  the  error  conmitted  by  Htopping  after  any  term  Is  lees 
In  magnitude  than  the  first  term  neglected. 

Evaluation  of  Pr/~ &>  x_7  for  any  m. 


oo  m-1  oo 


e-2yty2  2 dy  dt< 


J 

* (1+t)  2 , then 


Let  H(t) 


H(t)  . x 

V*>  ■ 7: 

H(rtl)(7)  - (. 

Then  | H(r+1) 

I V‘>l  S 

Therefore 

\ Pr/Tl,m>x- 


♦ Rr(t),  where 


f (*-y)r  H^r+1^(y)  dy,  and  where 
0 


l)r+1  * 2 r(r-t|)  (1+t)  * 2\  Then 
( J)\  < U(r+1)(0)|,  eo  that 


tr+1 

(r+1) ; 


B(r+1)(0)|  . 


J-o  J>  f+j 

zc 


I ,(x)  I 

J+l 


r(r<j)r(§«r+l) 

m 

_ rt+r+l 

nr(|)(r+l)i22 


I , (x),  where 

ri 


58 


J*0 , 1,  « • i 


Hence,  the  error  committed  is  leee  In  magnitude  than  the  first 
term  neglected.  We  find  that 

-x  -(*|> 

1 iW-5 — f - A I i(x)  >0,1,...  . 

H H H H 

By  repeating  the  above  recurrence  relation  several  times,  ve 

oo 

find  that  all  I 1(x)  can  be  made  to  fall  back  on  J e^t1/2  dt, 
3+2  / 

which  ve  can  get  from  a table  of  the  incomplete  gamma  function. 

Briefly  summarizing,  then,  ve  have  shown  hov  to  evaluate 

any  gfl  ^(t)  and  Pr/*^  m>  x_7  so  that  ve  can  use  3,1.6  and  3.1.7 

to  find  any  c.d.f. 

Inequalities. 

In  vhat  follows,  under  this  heading,  ve  shall  discuss 
certain  inequalities  related  to  the  distribution  of  the  difference 
between  tvo  Independent  chi-squares  considered  in  the  preceding 
sub-section,  3.1,  of  section  3. 


1 


oo 


Jt  ».! 

ft  ft  1 


“■■'•’■Ms  / 


JL  J 

and  since  a < (1-*^)  2 < 1,  it  follows  that 


JJ  jn  a 

3.1.8  2 2 h1(t)(l4^)  2<  g^ffl(t)  * 2 2 hx(t) 


and 


3.1.9  2 


£ 1-H1(x)__7  > Pr^”Ti>ni^x-7i  2^(1+ 


x 

whara  H^x)  . J hn(t)  dt, 

0 


x >0. 


If  n,a>  2.  than  ainoa 


n-ta 


(y+t) 


m-2  5-1  §-i 

> (y+t)  y2  * y 2 


*y> 


. t>0, 

Jtt 

V l-Hj^fx)^, 


, It  follows  that 
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5.1.10 


n-HB 


.00 

J • 


r(f)r(f) 


2t 


•mtrt 


n-HB 


- 1) 


n-t®  . 

^ r(|)r(|) 


n . (n-2)y 

_ o"l  * ,5f“ 

Again,  since  1^  (l+£)  5 © , it  follows  frost  5. 1.1  that, 

for  t > (^jp) , 


ass 

3.1.11  2*  ^(t)  * g^a(t)  * 2*^  hn(t)  /l  - 


Now  lotting  t -■— » 00  in  3.1.6  and  3.1.11  we  have 
a 

3.1.12  «n, *(*>/*  2_5  *>n(t),  »o  that 

8q  a(t)  has  the  saae  order  of  contaot  at  +00  as  the  p.d.f.  of 
2 

X with  n degrees  of  freedoa.  Similarly,  gn  #(t)  has  the  saae 

2 * 

order  of  contact  at  -00  as  the  p.d.f.  of  X with  a degrees  of 
freedoa. 

In  3.1.10  letting  t~— » 0,  we  have 


in 


r(3=  - l) 

3.i.i3  %a(0)  ■ -gs- , and  .o 

2 2 " r(|)r(f) 


the  frequency  curve  does  not  meet  the  origin. 

The  oase  n«m. 

Under  this  hesdlng  ve  shall  consider  certain  properties 

of  the  p.d.f.  and  c.d.f.  of  T . In  this  case  the  p.d.f.  of 

n,n 

Tr.  n Is  the  same  as  that  of  the  sum  of  products  of  pairs  of 

Independent  N(0,  1)  variates. 

On  the  p.d.f.  of  T 
SxS 

The  p.d.f.  Is  symmetric  and  so  we  shall  consider  g^  n(t) 
for  t > 0.  If  we  put  n«*m  in  3.1.13,  23  get 


3.1.1^ 


r(n-l) 


r(ffi 
2r(|)r(|)  ' 


We  can  show  that  gj^  n(0)  is  a decreasing  function  of  n. 
Differentiating  3,1.14  ve  have 


2r(|)r(|)  r'(|) 

-r-  fiLl  (0)  r — — 

r(n^l)  %n  r(»gi)  r(|) 
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Frcei  Cramer  (£ 5__7>  P«  131) , we  have  that  r'(n)/r(n)  ie  an 
increasing  function  and  so  gn  n(0)  is  a decreasing  funotion  of  n. 
That  Is, 


3.1.15 


k «n,n«»  < °* 


Finally,  it  is  easy  to  show  that 

3«1. l6  ^ g^  n(t)  ■ 0 at  t=0,  and 


consequently  that  ^ n(t)  has  its  maximum  at  the  origin. 


If  in  3.1.1  vs  put  n*m,  end  let  2x  = t(y-l),  we  get 


3.1.17 


If  n*l,  g^  ^(t)  « jj  Kq( t).  It  la  known  that  KQ(t)  ie  asymptotic 

to  both  axes  and  has  a logarithmic  singularity  at  the  origin. 
Using  a well-known  expansion  for  KQ(t),  £ lh£/ , we  find  that 

«iog!  • 


3.1.18 


as  t 


0. 


The  moment  generating  function  of  T _ becomes 

n,n 

n 

o * 3 

M(e)  - (1-0  ) and  so 


It  la  worth  noting  that  If  ve  expand  TQ  Q and  use  the  known 

p 

moments  of  X , we  gat  as  a by-product  that 


where  s*0,l,, . . and  n*l,2, ... 

On  the  c«d«f«  of  T » 
n,a 

For  large  values  of  n we  nay  wish  to  use  the  normal 
approximation,  since  T la  asymptotically  normal  (0,  */n).  It  la 

n 

here  appropriate  to  use  a result  due  to  Berry  /"  3 ].  Let  X*  L x, , 

1*1  1 

2 

where  the  Xj  are  Independent  r.r. *s.  Let  EX*  o , Var  X*a  , the 
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0(x) 


dt,  then 


-oo 


Sup 

-ocK  x<oo 


T(x)  - 0(2-2—)) < 


vhere 


(Sir) 


< C < 1.88,  according  to  Berry. 


If  we  let  » ■|{X^  - Y^),  then  we  may  write  Tn  n * s 

where  X^  and  are  independent  N(0,  1)  variates.  Putting  n-m=l, 

and  2x  ■ t(eec  © - 1)  in  3.1.1,  we  find  that  E | U^|  ^ » 8/n. 

If  F„(x)  Is  the  c.d.f.  of  T . then 
n n,n' 


3.1.19 


Sup  I Fn  (x/n)  - G(x) 
-oo<  x<  oo  ' 


8C 

n^n 


Consider  next 


oo 


x 


tn  ^(t) 


dt. 


Integrate  by  parts  where 


4 


^5 

U ■ t K^t),  dV  ■ tn“*  fit,  and  uae  the  relation  t K^(t)  ■ 
nK^(t)  - t K^ft),  getting 

oo  oo 

J tn+1  K^ft)  dt  - xn+1  Kjx)  + (2n+l)  J tn  Kn(t)  dt. 

x x 

Replacing  n by  (n-l)/2,  we  have 

ia-li 

o /x/2) 

3.1.20  pn<a(x)  - »„<*>  - V«A  (I) 

- Vx>  - f <v<x> 

-X  ,-x 

Hence,  knowing  that  gg  2(x)  " “2"  > S^j/x)  ■ ip  U+x), 

-x 

F2(x)  « 1 - -2^-  , we  can  obtain  ell  Fg^gU).  Similarly,  having 
a table  of  the  Beeeel  functions  Kq(x),  K^(x),  . . . we  can  get  all 

F2n+1^'  lf  we  taow  Fi^’ 

Evaluation  of  F^x). 

00 

We  have  %fl  - Fl(x)_7  ■ J ^(t)  dt.  Ueing  tables  of 


% 
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KQ(t)  ve  oould  evaluate  F^(x)  by  numerical  quadrature  j however, 
ve  give  en  alternative  method  here.  Integrate  by  parts  where 
U » l/t,  dV  ■ tKQ(t)  dt,  and  use  the  relation  « - 

-t^ft),  n«l,2,...,  getting 

oo  oo 

J dt  * Ei(l)  • / Ki<t)  r • 

X X 

If  we  carry  out  this  integration  by  parts  repeatedly,  we  find  by 
induction  that  we  get  an  alternating  series,  in  which  the  error 
committed  by  stopping  with  a given  term  is  less  in  magnitude 
than  the  first  term  neglected.  If  dy  is  the  r-th  term  of  the 
series,  then 


(-l)r+1  2r 

r (h 


r(r+|) 


r-l 


Kr(x), 


r*l,2, . • 


Furthermore,  if  s and  k are  any  two  positive  integers. 


2s- 1 2k 

Z d > it  jTl  • F- (x)  7>  Z d . 
1r-  i--1r 


3.1.21 


3.2  Special  case  2, 


In  this  section  vs  shall  give  an  expression  for  the  p.d.f. 
of  an  Indefinite  q.f , when  the  latent  roots  of  the  matrix  of  the 
q.f.  are  equal  in  pairs. 

Theorem  3.2. 

Let 


+ \ Yk  “ ak+l  Yk+1 


...  - an  Yn)> 


vhere  the  ai>  0,  end  the  Y^  are  Independent  r.v. 'e  each  having 

a chi-square  distribution  with  two  degrees  of  freedom.  If  f(q) 
le  the  p.d.f.  of  then 


k 

f (q)  » Z e 
J-l 


1- 


n-2  k 

a.  TT  (a.-a.) 

J 1.1  J 1 


-1 


n 


TT  («,+».) 

e»k+l  J 


q>  0 


n a.  n-2  k , n . 

- Z e 3 a.  TTCa  -fa.)"1  TT  (a«-a  ) , q<0  . 

J-Jc+1  0 i-1  J 1 s-k+1  J B 


Proof. 


The  moment  generating  function  of  ^ is 


M(t)  * £ (l-a^t) . . . (l-a^t)  (1+a^t) . . . (l4«nt)__7  -1.  Then 


oo 

f (q)-^j  f e-it<1  M(it)  dt. 
-oo 


We  can  evaluate  f(q)  using  contour  integration.  Let  us  take  as 
our  contour  in  the  complex  plane,  the  real  line  from  -R  to  +R 
and  then  the  semi-circle  of  radius  R,  in  the  lover  half-plane  if 
q>  0,  and  in  the  upper  half -plane  when  q< 0.  In  both  oases  the 
value  of  the  integral  around  the  semi- circle  tends  to  zero  as 
R — » oo,  if  n>  2.  Hence 

f (q)=2iti  £ Sun  of  residues  of  the  integrand  at  t®-—-,  J=l,...,k  7 

J 

q>  0,  and 

■2iti  /"Sum  of  residues  of  the  integrand  at  t=~,  J*?k+1, , . .,n_7 

J 

0. 


Evaluating  the  residues,  ve  get  the  form  stated  in  the  theorem. 


CHAPTER  IV 


FURTHER  BOUNDS  ON  THE  C.D.F.  OF  A DEFINITE  QUADRATIC  FORM 
IN  INDEPENDENT  IDENTICALLY  DISTRIBUTED  VARIATES  - EACH 
CENTRAL  NORMAL  OR  MORE  GENERAL 

4.1  The  central  normal  case. 

If  ve  make  uee  of  the  p.d.f.  of  Q^  when  the  latent  roots 
of  the  matrix  of  Q^  are  equal  in  pairs,  ve  can  obtain  some 
convenient  bounds  on  the  p.d.f.  of  Q-  when  the  latent  roots 
are  not  necessarily  equal  in  pairs.  Let 

Q2m  3 i<aixl+*  * •+%4+am+l4+l+*  * •+*2a4n)'  whero  ai>  °'  *■!*..  «i2»i 

end  where  , o2=*am+2,  . . . , am*a2m»  tlxen  th®  moment  generating 

function  of  is 

M(t)  ■ /*l-a1t)...(l-ajnt)-7  and  hence  the  p.d.f. 

oo 

of  Q^  is  h(q)  » ^ f e’lt<l  M(it)  dt. 

-oo 


Using  the  calculus  of  residues  ve  readily  find  that 
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4.1.1 


. SL 

h(l)  - £ e J a“"2  TT  (a.-a J"1 
>1  0 k-1  J K 


Suppose  now  that  we  want  the  p.d.f.  of 


^2n  “ |(aixi+* 


,+a2nW‘ 


where  a^  > > 


> a^  >0. 


We  form  the  two  expressmens 


%-\£ «i(^)  + '3«3+^)  + • ••  - e2n-1^2n-l+  4a )J>  “* 


- I ♦ a kiZ*+i£>  ♦ ...  ♦ a2a(X^n.1+  X^)_7,  so  that 


, >2  ,t2 


S-  ^2n  - * 

We  can  find  the  p.d.f.  for  and  by  using  4.1  so  that  we  have 
bounds  on  the  c.d.f.  of  Q^.  Let  f^q),  f(q),  f^(q)  4e  the 
p.d.f. 's  of  Qy,  Qgn,  Qt  respectively.  Then 


r 

J tjflq)  dq  < Pr/Q^l  tj  < f t L (q)  dq,  where 
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*„(«) 


n ft 

£ o 
J«1 


8J_1  jj  - 

tfj 


. _a_ 


fT(q)  - £ o ***  a?;1  TT  (®o4  - a2k)“1 


J-l 


2j  k:iV  2J 

WJ 


Remarks, 

(I)  The  above  inequality  was  suggested  by  Professor 
Hotelling. 

(II)  The  above  method  could  be  extended  to  cover  the 
case  of  an  indefinite  q.f. 

k,2  The  general  case. 

In  this  section  we  shall  discuss  briefly  a system  of 
inequalities  for  the  distribution  of  a definite  q.f.  Let 

Si  * |(®i*i  * •••  * “bv*  yher“ 


°_7  * P|  i*l|  see  Then 

Pr(l?  <c,...,I^<c)  * (l-p)n.  It  follows  that 
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Pr  •••  + »nL7>  (1-P)°  • 


Let  a1  > > . . . afl  ^ 0,  getting  the  following  system  of 

inequalities: 


an-7  < (i'P)11  , 

Pr  i"®n<!(*a.i+  anL7  < (l-p)“  + (")(l-p)n‘lp, 

Pr  /«„<  |{»n.2t  en.lt  an)_7  < (X-p)n+(J)(l-p)I“V(”)(l-P)n' V, 


Pr  /«n<  |(tt1+...+aa)_7  < l-p" 


The  above  system  was  suggested  to  the  writer  by  Professor 
Harold  Hotelling.  The  following  Improvements  are  due  to 
Professor  S.  N.  Roy. 

2 

Suppose  that  the  distribution  of  is  known  and  that 

*<*?  2 r>  ■ V 


2c 


a . +a 
n-l  n 


) • P 


2* 


2c 


an-2+an-l 


I 


2c 


l n 


n 


Then  It  follows  that: 

(1_Pn)n<  P*^  <c)  * (1-p^11, 

(l-pn)n<  Pr^o)  S (l-p2)n  + (“)(l-p2)n_1  p2, 

(i-Pn)n<  d-P5)n+  (J)(i-P3)n'lP5+  (?)d-P3)n-2  P3, 

(1-Pn)n<  Pr(Qn<o)  < l-p“  . 

Hence  we  have  one  lower  bound  and  a whole  system  of  upper 
bounds.  Obviously,  we  would  want  the  least  upper  bound,  in 
practice.  Incidentally,  P^<  Pg<  • • • < pQ. 
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